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(5) INTRODUCTION 

The thrust of the proposed, “Breast Cancer Screening Using Photonic Technology” research 
project is to develop a breast cancer screening modality that uses noninvasive and non-ionizing 
near-infrared (NIR) light for imaging and diagnosis of cancerous lesions of human breast. The 
imaging method involves illuminating the specimen with ultrashort NIR pulses of laser light and 
construction of images using two approaches. The first, known as the shadowgram method, 
utilizes the image bearing component of the forward-transmitted light to form direct shadow 
images. The second, referred to as the three dimensional (3-D) inverse reconstruction 
method, makes use of the measured transmitted and scattered-light intensity profiles, known or 
estimated optical properties of the sample, a model for light propagation through turbid media 
and a sophisticated computer algorithm to construct images of the interior structure of the 
specimen. 

We use a time-sliced detection technique that records a time sequence of two-dimensional (2- 
D) spatial intensity profiles of light transmitted through the specimen. These 2-D spatial profiles 
provide direct shadow images for thin samples. Alternatively, and for thick samples, these 2-D 
profiles are used to reconstruct 3-D images that provide information about the location of the 
object as well. Tuning the wavelength of imaging light to a known absorption resonance of 
adipose tissue in a specimen enabled highlighting of the adipose tissue from the fibrous tissue in 
the specimen. This result demonstrates the diagnostic potential of the spectroscopic imaging 
method. 

(6) BODY 

The research during the current reporting period proceeded along both the parallel 
approaches, i.e., the shadowgram method, and the inverse reconstruction technique, and we 
have major achievements to report in both the arenas. The tasks performed and the progress 
made may be broadly grouped as follows: 

6.1 Setting up and testing of experimental arrangements, 

6.2 Development and implementation of the inverse image reconstruction method, 

6.3 Time-sliced shadowgram imaging with excised human breast tissues, 

6.4 Spectroscopic imaging of excised human breast tissues, and 
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6.5 Demonstration of 3-D inverse image reconstruction. 


We will briefly outline our accomplishments in each of the above five areas, and refer to 
appended publications (Appendices 1-7) for detailed description where applicable. 

6.1 Setting up and Testing of Experimental Arrangements 

The major thrust of research in the first year of the approved Statement of Work (SOW) 
was to set up the experimental arrangements for time-sliced imaging [Technical Objectives 
(TO) 1, 2 and 3; Tasks 1-5]. We have set up the following experimental arrangements for 
near-infrared shadowgram imaging. 

(a) Time-sliced Imaging Arrangement with an Ultrafast Electronic Time Gate : As 
outlined in TO 3 (Tasks 5, and 6) of SOW, we have completed the setting up and testing of 
a time-sliced imaging arrangement. What is even more important is that, we have used the 
arrangement to obtain time-sliced images of normal and cancerous human breast tissues, 1 ’ 2 
as well as to reconstruct 3-D inverse images of objects in highly scattering medium. 3,4 

The arrangement uses the output (750-850 nm) of a femtosecond Ti:sapphire laser and 
amplifier system 3 generating 150 fs, 1 kHz repetition rate, 5 mJ pulses for sample 
illumination. Time-sliced 2-D images are recorded by an ultrafast gated intensified camera 
system (UGICS). The UGICS is a compact gated image intensifier unit that is fiber- 
optically coupled to a CCD camera. It provides an electronic time gate whose FWHM 
duration may be adjusted to a minimum of approximately 80 ps and the gate position may 
be varied over a 20-ns range in a minimum step size of 25 ps. Details of the arrangement 
is presented in Section 2 of Reference 3 which is included as Appendix 3 of this report. 

(b) Time-sliced Imaging Arrangement with Optical Kerr Gates : Optical Kerr gates 
(OKGs) hold the promise for providing temporal resolution of a few hundred femtoseconds, 
the limit being set by the pump pulse duration or the response time of the Kerr medium, 
whichever is the longest. In order to exploit this high temporal resolution we proposed 
development of combined OKG, polarization and Fourier space-gated imaging systems (TO 
1, Tasks 1-3). We have set up and tested (TO 2, Task 4) two OKG-based imaging systems. 
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The first system uses the second harmonic output (400 nm in use but tunable over 375- 
425 nm) of the Ti:sapphire laser system mentioned above to open the Kerr gate and the 
fundamental output (800 nm in use, but tunable over 750-850 nm) for imaging, and CS 2 as 
the Kerr medium. A schematic diagram of the system is provided in Appendix 5. 

The second system uses the second harmonic output (640 nm in use but tunable over 
610-650 nm) of our proprietary Cnforsterite laser and regenerative amplifier system 
mentioned above to open the Kerr gate and the fundamental output (1280 nm in use, but 
tunable over 1220-1300 nm) for imaging, and CS 2 as the Kerr medium. The laser system 
provides 200 fs, 1 kHz repetition-rate, 50 pJ pulses. A schematic diagram of the system is 
provided in Appendix 6. 

Our comparison of the OKG and the UGICS based systems indicates that while the 
OKG-based system provides higher temporal resolution, the UGICS-based system is more 
user friendly, easier to implement and provides higher signal since integration over a wider 
time window is employed. 

(c) Spectroscopic Imaging Arrangement : We have assembled a spectroscopic imaging 
arrangement to explore the diagnostic potential of optical biomedical imaging techniques. It 
uses the 1210-1300 nm continuous-wave mode-locked output of a Cr.forsterite laser for 
sample illumination. A Fourier space gate 6 and a polarization gate 7 are used to select the 
image-bearing photons and discriminate against the multiple-scattered photons. The 
shadowgram images are recorded by a 128x128 pixels InGaAs NIR area camera (Sensors 
Unlimited SU-128). The arrangement is detailed in Section 2 of Reference 2 (Appendix 
2) and a schematic diagram is presented in Appendix 7. 

6.2 Development and Implementation of an Inverse Image Reconstruction Method 

We have developed 3 a novel inverse image reconstruction (HR) method that makes use 
of the sequence of time-sliced 2-D transmitted light intensity distribution measured by the 
experimental arrangements detailed in Section 6.1 to reconstruct 3-D images of object inside 
highly-scattering turbid media (TO 4, Tasks 8, 9). Diffusion approximation of the radiative 
transport theory 8 is used to model the light propagation through the turbid media. The 
formulation of the forward problem assumes a slab geometry described in cylindrical 


7 



coordinates. A Green’s function perturbative approach under the Rytov approximation 9 is 
used to produce a linear inverse algorithm. The use in the inversion algorithm of a 2-D 
matrix inversion with a 1-D Fourier-transform inversion based on the cylindrical symmetry 
greatly reduces the computation time. The algorithm has been implemented on a Silicon 
Graphics Origin-2000 computer equipped with 4 parallel 195-MHZ CPU’s. The inverse 
reconstruction method is detailed in Appendix 3. 

6.3 Time-sliced Shadowgram Imaging with Excised Human Breast Tissues 

The efficacy of the time-sliced imaging technique outlined in Section 6.1 was used to 
obtain shadowgram images of excised human breast tissue specimens (TO 5-7, Tasks 12, 
14). The specimens were obtained from our collaborators at the Memorial Sloan Kettering 
Cancer Center and National Disease Research Interchange under an IRB approval at the City 
College of New York. Specimens with normal and cancerous tissues were used in the 
experiments. Specimens with infiltrating ductal carcinoma, infiltrating lobular carcinoma, 
and mucinous carcinoma were investigated. 

One of the key results of the study was that light transited through the cancerous tissues 
faster than through the normal tissue. Consequently, images obtained with light in the earlier 
time slices highlighted the cancerous tissues, while those obtained with light in the later time 
slices highlighted the normal tissues. Time-sliced imaging can thus separate out normal and 
cancerous tissues in excised specimens. Details of the experimental procedure, results, 
and time-sliced images are presented in References 1 and 2 (Appendices 1 and 2). 

6.4 Spectroscopic Imaging of Excised Human Breast Tissues 

In order to explore the diagnostic potential of optical imaging modalities, (TO 6, Task 
13) we carried out spectroscopic imaging of excised human breast tissues using the 
spectroscopic imaging arrangement mentioned in Section 6.1(c). We used excised breast 
tissue specimens with adipose and fibrous tissues. Shadowgram images were recorded by 
tuning the Cr:forsterite laser output to (a) 1225 nm close to an adipose tissue absorption 
resonance, 10 and (b) to wavelengths away from the resonance, such as, 1235, 1245, 1255, 
1265, 1275, 1285 and 1300 nm. Marked enhancement in contrast between the adipose and 
fibrous tissues was observed for the near-resonant image recorded using the 1225-nm light. 
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The experimental procedure, results, and spectroscopic images are detailed in 
References 1 and 2 (Appendices 1 and 2). Results of this spectroscopic imaging 
experiments demonstrate the diagnostic potential of optical imaging. 

6.5 Demonstration of 3-D inverse image reconstruction 

The ability of the inverse reconstruction method outlined in Section 6.3 to generate 3-D 
images of an object embedded in a turbid medium using the 2-D time-sliced transmission 
measurements outlined in Section 6.1(a) was demonstrated (TO 4, Tasks 8, 9). The 
absorbing object was a 3x3x10 mm 3 rectangular parallelepiped made of aluminum and 
painted black. The scattering medium was Intralipid-10% in water with a reduced scattering 
coefficient, p s ’ of 0.4 mm' 1 and an absorption coefficient p a of 0.02 mm' 1 . The Intralipid 
suspension was contained in a 60 mm long (approximate size of an averaged compressed 
human breast), 200 mm diameter cylindrical plexiglas cell. The object was suspended on the 
axis at a distance of 15 mm from the exit surface of the cylindrical cell. Our goal was to 
obtain 3-D images of the embedded object from time-sliced transmitted intensity 
measurements. 

A sequence of 2-D slices at 100 ps intervals were recorded over a 5-ns range. Three 
dimensional images of the object with good lateral spatial resolution were obtained. Spatial 
resolution in the axial direction needed improvement. Simulation using data with different 
levels of Gaussian distributed noise revealed that noise level in the experimental data may 
have contributed to the spreading out of the reconstructed image along the axial direction. 
We are working on reducing the noise level in the experiment and to make the reconstruction 
algorithm more noise resistant. The inversion technique and the reconstructed images are 
detailed in References 3 and 4 (Appendices 3 and 4). 

(7) KEY RESEARCH ACCOMPLISHMENTS 

• Demonstrated that time-sliced near-infrared optical imaging enables highlighting of 
normal and cancerous tissues in excised breast specimens. 

• Near-infrared spectroscopic imaging can highlight different types of tissues (such as, 
adipose and fibrous) when the wavelength of light is tuned to an optical resonance of 
any of the tissues. 
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• Results of near-infrared spectroscopic imaging on excised tissues illustrate the 
diagnostic potential of optical imaging. 

• Reconstruction of 3-D images of objects in turbid media from 2-D time-sliced 
transmission was demonstrated. The method has the potential to provide location of 
the object in 3-D that is not available from normal x-ray mammography. 

• Spectroscopic and time-sliced measurements together have the potential to provide 
more diagnostic information than other techniques, such as, x-rays. 

(8) REPORTABLE OUTCOMES 
Articles 

1. S. K. Gayen and R. R. Alfano, “Sensing lesions in tissues with light,” Opt. Express 4, 
475 (1999). 

2. W. Cai, S. K. Gayen, M. Xu, M. Zevallos, M. Alrubaiee, M. Lax, and R. R. Alfano, 
“Optical tomographic image reconstruction from ultrafast time-sliced transmission 
measurements,” Appl. Opt. 38,4237 (1999). 

3. S. K. Gayen, M. E. Zevallos, B. B. Das, and R. R. Alfano, “Spectroscopic and time- 
sliced imaging of human breast tissues,” SPIE 3597, 508 (1999). Proceeding of the 
SPIE Conference on Optical Tomography and Spectroscopy of Tissue III, 24-28 January 
1999, San Jose, California. 

4. M. Xu, S. K. Gayen, W. Cai, M. E. Zevallos, M. Lax, and R. R. Alfano, “Time-sliced 
three dimensional inverse image reconstruction of objects in highly scattering media,” 
SPIE 3597, 2 (1999). Proceeding of the SPIE Conference on Optical Tomography and 
Spectroscopy of Tissue III, 24-28 January 1999, San Jose, California. 
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Tomography and Spectroscopy of Tissue III, 24-28 January 1999, San Jose, California. 

2. S. K. Gayen, M. E. Zevallos, B. B. Das, and R. R. Alfano, "Near-infrared spectroscopic 
and time-sliced imaging of human breast tissues." Presented at SPIE's International 
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Symposium on BiOS '99, Conference 3597: Optical Tomography and Spectroscopy of 
Tissue III, 24-28 January 1999, San Jose, California. 

3. S. K. Gayen, M. E. Zevallos, and R. R. Alfano, "Time-resolved transillumination 
imaging of normal and cancerous human breast tissues with a picosecond electronic 
time gate, " Bull. Am. Phys. Soc. 44, 117 (1999). Paper BC32 7 presented at the 
Centennial Meeting of the American Physical Society, 20-26 March 1999, Atlanta, 
Georgia. 

4. S. K. Gayen, M. E. Zevallos, and R. R. Alfano, "Time-sliced and spectroscopic near- 
infrared imaging for cancer detection." Paper CFL4 presented at the Conference on 
Lasers and Electro-Optics (CLEO ’99), May 23-28, Baltimore, Maryland. 

Degrees Obtained and Employment Received 

The Ph. D. research of Manuel E. Zevallos (U.S. Citizen, Hispanic) entitled, “Near 
Infrared Optical Imaging and Light Propagation in Highly Scattering Random Media” 
(defended February 11, 1999) was supported in part by this award and the NASA 
Institutional Research Award Program. Dr. Zevallos is now employed with IBM. 

(9) CONCLUSIONS 

The work carried out just in the first year of the award leads to some useful conclusions. 
First, time-sliced imaging measurements reveal that light transits through the cancerous 
tissues faster than the normal tissues. This feature enables highlighting of cancerous and 
normal regions of tissues. Second, results of spectroscopic imaging measurements led to 
the conclusion that if any tissue had a specific spectroscopic signature (such as, enhanced 
absorption), then resonant images that exploit that signature can highlight that tissue from 
other types of tissues. Finally, time-sliced 2-D transmission spectroscopic measurements 
provide a wealth of information for reconstruction of 3-D tomographic images. 

“So What Section” 

The implication of the first conclusion is that time-sliced imaging offers the possibility of 
highlighting cancerous lesions in human breast. If the lesion is not too deep into the breast, 
a direct imaging using backscattering geometry may be possible. 

The second conclusion points to the most important potential of optical imaging, that of 
being able to diagnose the lesion while detecting it. X-ray mammography, the current 
screening method, can not diagnose cancer. 
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The final conclusion implies that using time-sliced 2-D transmission measurements 3-D 
tomographic images of human breast may be obtained. 

The research carried out so far lends credence to the original hypothesis that photonic 
technology has the potential to emerge as a noninvasive modality with diagnostic potential 
for breast cancer monitoring. Additional funding (federal and corporate) is needed to 
advance this novel technology. 
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Appendix 1 


Sensing lesions in tissues with light 


S. K. Gayen and R. R. Alfano 

Institute for Ultrafast Spectroscopy and Lasers, New York State Center for Advanced Technology for Ultrafast 
Photonic Materials and Applications, Physics Department, The City College and Graduate School of the City 
University of New York, New York, NY 10031 

gayen QPscisun.sci.ccny. cuny. edu, alfano QPscisun. sci. t en warn v.edu 


Abstract: The use of light for probing and imaging of biomedical media 
offers the promise for development of safe, noninvasive, and inexpensive 
clinical imaging modalities with diagnostic ability. Various properties of 
light together with the ways it interacts with biological tissues may provide 
‘multiple windows’ to peer inside body organs. Principles and methods for 
extraction of information about body functions and lesions that capitalize on 
temporal, spectral, polarization, and spatial characteristics of transmitted 
light are briefly outlined. As illustrations of the potential and efficacy of 
light-based techniques, time-sliced and spectroscopic images of normal and 
cancerous human breast tissues recorded with a femtosecond Tiisapphire 
laser and a broadly tunable Cnforsterite laser, respectively, are presented. 

©1999 Optical Society of America 

OCIS codes: (170.3880) medical and biological imaging; (170.6920) time-resolved imaging; 

(290.7050) scattering, turbid media; (999.9999) optical mamm ography 
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1* Introduction 
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^ PSI ba f S Wltnessed a ra Pid growth of research activities to exploit some of 
those attributes of light for the development of optical biomedical imaging and biopsy 

EST ' AW “ U * h “ ™Pr es sive array of body-imaging techniques, fuch as x-ray 

r“f 8 g ; y com P uted tomography, magnetic resonance imaging, ultrasound and 

radioisotope imagtng are currently available to yield useful information [2] there are 
important limitations of safety, resolution, cost, and lack or limited specificity to key 
chemicals necessary for functional body monitoring. More importantly, convention2 
hmques sometimes fail to provide necessary information. For example x-rav 
mammography, the current ‘gold standard’ for monitoring breast cancer, is not suitable S 

O^Sl 8 im enSe ’ ’Tf r eaStS ’ “ ,d may DOt distklguish ^tween malignant and benign mmors 
Optical imaging modalities are being pursued to alleviate most of these limitations and to 

' ^agnostic ability through the exploitation of the spectroscopic characteristics of 
ssues. Consequently, a number of continuous-wave, time-resolved, and frequency-domain 

X of riwf T S 5 aVe CVOlVed ° VCr * e yearS [3 - 16 ’^' lt should afro be noted that 
any of the above-mentioned imaging techniques provide one type of image only, while bv 

c angrng the wavelength of light one can obtain multiple images that highlight the various 
molecular components. In addition, light interacts with biological tissues in a number of 
different ways providing different means to probe the body. 

In this article, we first present a brief overview of the various light-tissue interactions and 
the properties of light that enable sensing the interior of body organs, such as breast and 
prostate. We then illustrate the biomedical sensing potential of light through the time-sliced 
optical imaging and spectroscopic imaging of normal and cancerous human breast tissues. 

2. Background 

Shadowgram optical imaging was first attempted in 1929 by Cutler [3] who used white light 
lluminate die breast and look for pathology in the transmitted light. Inadequate spatial 
resolution attainable with the available light sources of that time, and the rapid Sogress in x- 
ray mammography eclipsed the development of light-based method, that was afroltnown as 

rSr," Cartography [2]. With the advent of broadly tunable nL2Z3 
(NIR) solid-state lasers, such as Ti:sapphire and Cr: forsterite, and the development of charge- 
coupled device (CCD) and NIR area cameras, there has been a surge in research activities for 
the development optical body-scanning modalities since the mid-1980’s. 
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The purpose of an optical body-scanning modality is to identify, locate, and diagnose a 
‘lesion’, such as a tumor, inside human body. The procedure involves illuminating the 
relevant part of the body, such as the breast or the prostate, with bright light of appropriate 
wavelength and search for indication of pathology in the emergent light. The physical basis 
for the development of an optical modality is the difference between the interaction of light 
with the lesion and the surrounding tissues. Major problems that impede the development 
include attenuation of light due to scattering and absorption by the tissues, and image blurring 
due to scattering. The nature and extent of the changes in the properties of light that result 
from its interaction with biological tissues may be different for normal and diseased tissues, as 
well as for different types of tissues, and thus provide bases for the development of optical 
imaging and optical biopsy methods. 

The high spectral brightness of laser light makes it the most common choice of researchers 
for optical biomedical imaging. The enabling attributes of laser light include wavelength, 
coherence, polarization, pulse duration in the case of mode-locked and Q-switched lasers, 
intensity and directionality. Each of these properties may change as a result of light-tissue 
interaction that depends on the wavelength of light. A proper choice of wavelength(s) is a 
crucial consideration for developing an optical imaging modality. Light in the ‘therapeutic 
window’ of 700-1300 nm is less absorbed by the tissue constituents compared to near¬ 
ultraviolet, visible, or infrared light and is commonly used for tissue imaging. 

The common linear interactions of a beam of light with biological tissues include 
scattering, absorption, specular and diffuse reflection, and diffuse transmission. Scattering 
changes the directionality, intensity, coherence, polarization, and pulse-width, of the incident 
light. Consequently, light transmitted through a scattering medium breaks up into ballistic , 
snake , and diffuse components [5,9,1]. The early-arriving, forward-propagating ballistic and 
snake photons carry image information while the multiply-scattered diffuse light blurs the 
image, and in extreme cases buries it in the background noise. The difference in transit time 
and changes in the above-mentioned characteristics are exploited for selecting the image- 
bearing light to form direct transillumination images, or, ‘shadowgrams’. Since the image¬ 
bearing photons change the least, the idea is to devise a gate that will let a significant fraction 
of the photons with a specific initial property through but block others. Realization of this 
idea has led to the development of time , coherence , polarization , nonlinear optical , and space 
gates [4-16,1]. However, for a highly scattering and thick medium the image-bearing light 
becomes too weak to form a shadow image. One then resorts to inverse methods that make 
use of the scattered light intensity patterns measured around the object, a mathematical model 
for description of light propagation in scattering media, and a sophisticated computer 
algorithm to reconstruct an image of the object [17]. 

Fluorescence that results from the absorption of light by a molecular constituent, such as 
tryptophan, tyrosine, collagen, elastin and flavin, of a tissue is a characteristic of that 
constituent. Changes in the fluorescence spectra, lifetime, intensity, and intensity ratios for 
different excitation wavelengths may provide markers for lesions and facilitate their imaging 
[18]. Incident light may be frequency shifted due to Raman scattering by the vibrational 
modes of the molecular constituents of a tissue and differences in the Raman spectra of 
normal and diseased tissue may provide diagnostic ability [19]. Nonlinear optical interactions 
of two photon absorption-induced fluorescence and second harmonic generation have also 
been used to generate histopathological maps of biological tissues [15,16]. 

The organs of human body with potential to be investigated by optical spectroscopic 
techniques include breast, prostate, brain, bladder, bone, cervix, colon, eye, digestive and 
gynecological tracts, skin, and teeth [20]. The requirements of resolution, penetration depth, 
diagnostic ability, and specificity to key chemicals of an optical imaging modality depend on 
the organ and the lesion being investigated. 

The remainder of this article focuses on our NIR time-sliced and spectroscopic imaging of 
human breast tissues in vitro. Optical imaging of breast is complicated by the extremely low 
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ratio of image-bearing photons to multiply scattered image-blurring photons, and paucity of 
high contrast targets that are indicative of the presence of cancer and suitable for imaging 
Our thrust is to develop a better understanding of light transport through different types of 
tissues, obtain NIR transillumination images, and explore the feasibility of distinguishing 
between normal and cancerous tissues from those images. 5 

3. Materials and Methods 

We have mvestigated excised female human breast tissue specimens with normal fibrous and 
a lpose tissues, as well as normal and cancerous tissues. Types of cancer included 
infiltrating ductal carcinoma, and infiltrating lobular carcinoma. The tissue specimens were 
placed between two glass plates and compressed to provide a uniform thickness The 
specimens were provided to us by National Disease Research Interchange and Memorial 
Moan Kettering Cancer Center under an IRB approval from the City College of New York 

wirtf W ^ P resent representa ^ ve results of a time-sliced imaging measurement on a specimen 
wiA infiltrating ductal carcinoma, and a spectroscopic imaging measurement on a specimen 
withachpose and fibrous tissues to demonstrate the potential of optical imaging methods. 

Hie time-sliced imaging approach is an extension of the idea of time-resolved early-light 
imaging except that picosecond-scale slices of die transmitted light pulse are used to record a 
sequence of two-dimensional (2D) images for different positions of the time gate The 
experimental arrangement made use of 800-nm, 130-fs duration, 1 kHz repetition-rate pulses 
from a Ti.sapphire laser and amplifier system for sample illumination, and an ultrafast gated 
in ensified camera system (UGICS) for recording 2D images. The UGICS comprised a 
compact time-gated image intensifier unit fiber-optically coupled to a CCD camera It 
provided an approximately 80 ps-duration gate whose position with respect to an external 
Q-igger pulse could be varied in a minimum step size of 25-ps over a 20 ns range. The average 
beam power used in the experiment was approximately 200 mW. The beam was expanded 
“ a , 3 -cmdmmeter central part of it was selected out using an aperture to illuminate the 
sample. The tune-sliced images were recorded by the CCD camera for different gate 
positions and displayed on a personal computer. 

The experimental arrangement for NIR spectroscopic imaging, detailed elsewhere r211 
made use of 1210-1300 nm continuous-wave mode-locked oJpufof a CnforS to to 

fl“ te . the Sa ™P Ie ; A F° uner s P ace gate [13] and a polarization gate [11] were used to 
select the image-bearing photons. A 50 mm focal-length camera lens placed on the optical 
axis at a distance of 50 mm from the aperture in the Fourier gate collected and collimated the 
w-spatial-frequency light filtered by the aperture and directed it to the 128x128 pixels 
sensing element of an InGaAs NIR area camera. Appropriate neutral density filters were used 
o main tarn the average optical power of the incident beam at approximately 35 mW for all 

t SSIST ima8 “ 8 “ Perim “ 1 ' ^ b ““ ™ P" *»W 

4. Results 

Time-sliced images of a 5-mm thick breast tissue sample comprising a piece of normal tissue 
and a cancerous piece with infiltrating ductal carcinoma for two different gate positions are 
displayed in the left frames of Fig. 1. The zero position was taken to be the time tfaTvdtf 

tltt P ro g a /' mm thiGk ghss cel1 fiIled with water. Corresponding spatial 
mtensity profiles integrated over the same horizontal area for all the images are shown in the 

n§ht ‘ i- ? e SaUent featUre ° f * e mages “ d “tensity profiles is that at early 
times (25 ps) more hght is transmitted through the cancerous region than the normal region 
whde at a much later time (275 ps) relatively more light arrives through the normal re^’ 
Lower scattering or/and higher absorption of light by cancerous tissues may account for the 
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observed temporal behavior. Since there is no significant absorption of 800-nm light by the 
breast tissues, we attribute the observed time-dependent difference in light transmission to the 
higher scattering of light by the normal tissue than by the cancerous tissue. Although we did 
not carry out independent measurements of the light transport characteristics of the tissues, 
our results are consistent with the literature values [22] of reduced scattering coefficient of 1.3 
mm"* for normal glandular tissues, and 1 mm" 1 for cancerous breast tissue at 800 nm. Images 
recorded with different time slices of the transmitted light could highlight normal and 
cancerous regions of an excised breast tissue specimen. 


Sample 



25 ps 



275 ps 




Fig. 1. Time-sliced 2D images (left) of a 25mmx9mmx5mm breast tissue sample with normal 
and cancerous components for gate positions of 25 ps and 275 ps. Spatial intensity profiles 
integrated along the same horizontal area are shown in the respective frames to the right. A 
schematic diagram of the sample is shown in the top left frame. The zero position was taken to 
be the time of arrival of the light pulse through a 5-mm thick glass cell filled with water. 


The thrust of the spectroscopic imaging experiment was to test if a spectroscopic 
difference between different types of tissues in a specimen provides any distinguishable 
signature in the transillumination image. We obtained images of a normal breast tissue 
sample comprising adipose tissue in the middle and fibrous tissues in the two ends with light 
of wavelengths in the 1225-1300 nm range. Some of the wavelengths were near-resonant 
with the adipose optical absorption band around 1203 nm, [23] while others were far removed 
from that resonance. Measurements were made both at the room temperature of 20 °C and the 
normal body temperature of 37 °C. The salient features of the images, displayed in Fig. 2, 
are: (a) adipose tissue appears much darker (less light transmission) than the fibrous tissues 
for the near-resonant 1225-nm image as compared to the off-resonance 1285-nm image; and 
(b) transmission of light through the specimen is higher at 37 °C than at 20 °C. Adipose tissue 
region appeared as a deeper trough in the spatial intensity profile of the 1225-nm image 
compared to that for the 1285-nm image. For a more quantitative description of the observed 
behavior, we monitored the image contrast, C(A,T) = (1 F - 1 A )/(I F - /„), where I a (A,T) is the 
optimal intensity value at wavelength A and temperature T on the spatial profile of the image 
at the adipose tissue location, and Ifi A T) is the corresponding intensity at the immediate 
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fibrous tissue region. Values of contrast at 1225 nm are 0.59 and 0.46 for 20 °C and 37 °C, 
respectively, while the corresponding values for 1285 nm are 0.25 and 0.10. As the laser 
output was tuned away from the absorption peak, the contrast between the adipose and fibrous 
regions in the images decreased. These results demonstrate that an appreciable spectroscopic 
difference may significantly enhance the contrast between different types of breast tissues in 
the transillumination image. The ratio of and/or difference between images recorded with 
resonant (or near-resonant) light and non-resonant light may enhance the image contrast even 
further. 



20 °C 37 °C 

(a) 1225 nm 



Pixels 


F 
A 

F 

20 °C _ 

(b) 1285 nm 




Pixels 


Fig. 2. Spectroscopic 2D images (left frames) of a 28mmxl2mmxl0mm normal human breast 
tissue comprising fibrous (F) and adipose (A) components recorded with (a) 1225 nm and (b) 
1285 nm light. Spatial intensity profiles of the images integrated over the same vertical area 
are shown in the respective frames to the right. 


5. Discussion 


The results of time-sliced imaging experiment demonstrate that optical imaging using optimal 
time slices of transmitted light is a promising method for imaging biomedical media. While 
the early-arriving light highlighted the cancerous tissue, the later-arriving light accentuated 
the normal component. Different types of tissues scatter light differently, and for specific 
sample configurations the later-arriving light may be as revealing as the early light. Another 
novel application of the sequence of time-sliced 2D images that is obtained by this approach 
is in fast three-dimensional tomographic inverse image reconstruction [24]. The time-sliced 
imaging approach enables fast acquisition of data over a wide frequency range that is 
necessary for 3D reconstruction, as compared to frequency-domain schemes that commonly 
use several discrete frequencies. 

Results of the NIR spectroscopic imaging experiments demonstrate the diagnostic 
potential of optical imaging. Combined time-sliced and spectroscopic imaging method has 
the potential to provide much more information than even the x-ray techniques. 
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Spectroscopic and time-sliced imaging of human breast tissues 
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ABSTRACT 

Two-dimensional transillumination images of in vitro human breast tissue specimens were recorded as a function of the 
wavelength of near-infrared laser light, as well as, the position of picosecond-duration temporal slices of transmitted light 
pulse. Spectroscopic imaging arrangement made use of 1225 - 1300 nm light from a chromium-doped forsterite laser for 
sample illumination, a Fourier space gate and a polarization gate to sort out a fraction of the image-bearing photons, and an 
InGaAs area camera for recording two-dimensional images. Time-sliced imaging approach used 800-nm, approximately 
130-fs duration, 1 kHz repetition-rate pulses from a Ti:sapphire laser system to illuminate the sample, and a gated imaging 
system that provided a variable-position, ~80 ps-duration electronic gate to record time-sliced two-dimensional images. 
Excised female breast tissue specimens comprising normal and cancerous, as well as, adipose and fibrous tissues were 
imaged Marked enhancement of image contrast between adipose and fibrous tissues in the specimen was observed when 
the wavelength of imaging light was near-resonant with an absorption band of the human fatty tissue centered on 1203 nm. 
In time-sliced images, earlier time slices highlighted cancerous tissues, while the later slices accentuated the normal tissues 
in most of the specimens. 

Keywords: Optical mammography; near-infrared absorption spectroscopy of tissues; spectroscopic imaging; 
transillumination; time-sliced imaging. 

1. INTRODUCTION 

The article presents the results of spectroscopic and time-sliced transillumination imaging 1 ' 3 measurements on excised 
human female breast tissue specimens. The thrust of the spectroscopic imaging experiments is to examine if a 
spectroscopic difference may be used to distinguish between different types of tissues in a specimen. Time-sliced imaging 
is an extension of the time-resolved early-light imaging except that different temporal slices of the transmitted light 
intensity profile are used to form images. Images recorded with different temporal slices of the transmitted light is found 
to selectively highlight different types of tissues, such as, adipose and fibrous or normal and cancerous in excised breast 
tissue specimens. 


2 . METHODS AND MATERIALS 

The experimental arrangement for near-infrared (NIR) spectroscopic imaging made use of 1210-1300 nm continuous-wave 
mode-locked output of a Ci* + :Forsterite laser to illuminate the sample. A Fourier space gate 4 and a polarization gate 5 were 
used to select the image-bearing photons and discriminate against the multiple-scattered image-blurring photons. A 50 mm 
focal-length camera lens placed on the optical axis at a distance of 50 mm from the aperture in the Fourier gate collected 
and collimated the low-spatial-frequency light filtered by the aperture and directed it to the 128x128 pixels sensing 
element of an InGaAs NIR area camera. Appropriate neutral density filters were used to maintain the average optical 
power of the incident beam at approximately 35 mW for all the wavelengths used in the imaging experiment. The laser 
beam was linearly polarized along the horizontal direction. 

The experimental arrangement for time-sliced imaging made use of 800-nm, approximately 130-fs duration, 1 kHz 
repetition-rate pulses from a self-modelocked Tirsapphire laser and amplifier system for sample illumination, and an 
ultrafast electronic gated imaging system (UEGIS) for recording two-dimensional images using picosecond-duration slices 
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of the light transmitted through the sample. 6 The UEGIS comprised a compact time-gated image intensifyer unit fiber- 
optically coupled to a charge-coupled device (CCD) camera. It provided an approximately 80 ps-duration gate whose 
position could be varied in 25-ps steps over a 20-ns range. The average beam power used in the experiment was 
approximately 200 mW. The beam was expanded by a beam expander, and a 3-cm diameter central part of it was selected 
out using an aperture to illuminate the sample. The transillumination signal recorded by the system at a particular gate 
position was a convolution of the transmitted light pulse with the gate pulse centered on the gate position. The time-sliced 
image was recorded by the CCD camera and displayed on a personal computer. 

Excised breast tissue specimens were provided to us by National Disease Research Interchange and Memorial Sloan 
Kettering cancer center under an IRB approval from the City College of New York. Samples with normal fibrous and 
adipose tissues, as well as normal and cancerous tissues were studied. Investigated cancerous tissues included those with 
infiltrating ductal carcinoma, and infiltrating lobular carcinoma with cytologic atypia. The tissue specimens were placed 
between two glass plates and compressed in an attempt to provide a uniform thickness. 

3. RESULTS 

To test if a spectroscopic difference between different types of tissues in a specimen provides any distinguishable signature 
in the transillumination image, we obtained images of a normal sample comprising adipose tissue in the middle and fibrous 
tissues in the two ends using light of wavelengths that varied from being near-resonant with the adipose optical absorption 
band around 1203 nm, 7 to being removed from it. Figures 1(a) and 1(b) show 'near-resonant images' recorded using 1225- 
nm light, and 'nonresonant image' recorded using 1285-nm light, respectively. Measurements were made both at the room- 
temperature of 20 °C and the body temperature of 37 °C. The salient features of the images are: (a) adipose tissues appear 
much darker (less light transmission) than the fibrous tissues for the near-resonant 1225-nm image as compared to off- 
resonance 1285 nm image; and (b) transmission of light through the specimen is higher at 37 °C than at 20 °C. Adipose 
tissue region appeared as a deeper trough in the spatial intensity profile of the 1225-nm image compared to that for the 
1285-nm image. For a more quantitative description of the observed behavior, we monitored the image contrast, C(X , T) = 
(If -IA)!(If where I^(X,T) is the optimal intensity value at wavelength X and temperature T on the spatial profile of 

the image at the adipose tissue location, and Ip{X,7) is the corresponding intensity in the immediate fibrous tissue region. 
Values of contrast at 1225 nm are 0.59 and 0.46 for 20 °C and 37 °C, respectively, while the corresponding values for 1285 
nm are 0.25 and 0.10. As the laser output was tuned away from the absorption peak, the contrast between the adipose and 
fibrous regions in the images decreased. For near-resonant images contrast was higher at 20 °C than at 37 °C. These results 
clearly demonstrate that an appreciable spectroscopic difference may significantly enhance the contrast between different 
types of breast tissues in a transillumination image. The ratio of and/or the difference between images recorded using 
resonant (or near-resonant) and nonresonant light may enhance the image contrast even further. 
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Figure 1. Two-dimensional transillumination images of a 28 mmxl2 mmxlO mm human breast tissue sample comprising adipose and 
fibrous tissues recorded using light of wavelength (a) 1225 nm, and (b) 1285 nm. Spatial intensity profiles of the images integrated 
over the same vertical area are shown in the rightmost frames. 
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Time-sliced transillumination images of a 5 mm thick breast tissue sample comprising a piece of normal tissue and a 
cancerous piece with infiltrating ductal carcinoma for gate positions of 25 and 275 ps are displayed in Figs. 2(a) and 2(b), 
respectively. The zero position was taken to be the time of arrival of the light pulse through a 5-mm thick glass cell filled 
with water. The spatial intensity profiles integrated over an horizontal area with a vertical width of 20 pixels are shown in 
Fig. 2(c). The cancerous and normal regions are clearly distinct in both the 25-ps and the 275-ps images. The spatial 
intensity distribution of the 25-ps image, displayed by the thin line in Fig. 2(c), shows much higher intensity values in the 
cancerous region compared to that in the normal region indicating much higher light transmission through the cancerous 
region at early time. Interestingly, the situation is completely reversed in the spatial intensity profile of the 275-ps image 
indicated by thick line in Fig. 2(c). Lower scattering or/and higher absorption of light by cancerous tissues may account 
for the observed temporal behavior. Since there is no known absorption by breast tissue at 800 nm, we attribute the 
difference in the relative light transmission between the normal and the cancerous human breast tissues at the two different 
times to the higher scattering of light by the normal tissue. The scattering characteristics of the infiltrating ductal 
carcinoma sample presented here is different from the highly invasive lobular carcinoma sample with cytologic atypia that 
we reported in Ref. 6. 


Cancer Normal Cancer Norma! 



(a) 25 ps (b) 275 ps 



Figure 2. Time-sliced transillumination image of the breast tissue sample described in the text for gate delays of (a) 25, and (b) 275 ps. 
(c) Spatial profiles of the normalized and integrated intensity distribution along a horizontal area of 20 pixel vertical width of the 25 ps 
(thick line) and 275 ps (thin line). 

In summary, both the spectroscopic and time-sliced imaging methods show tissue selectivity. A combined spectroscopic 
and time-sliced imaging approach has the potential to provide more information even than the x-ray techniques. 
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Appendix 3 


Optical tomographic image reconstruction from 
ultrafast time-sliced transmission measurements 


W. Cai, S. K. Gayen, M. Xu, M. Zevallos, M. Alrubaiee, M. Lax, and R. R. Alfano 


Optical imaging and localization of objects inside a highly scattering medium, such as a tumor in the 
breast, is a challenging problem with many practical applications. Conventional imaging methods 
generally provide only two-dimensional (2-D) images of limited spatial resolution with little diagnostic 
ability. Here we present an inversion algorithm that uses time-resolved transillumination measure¬ 
ments in the form of a sequence of picosecond-duration intensity patterns of transmitted ultrashort light 
pulses to reconstruct three-dimensional (3-D) images of an absorbing object located inside a slab of a 
highly scattering medium. The experimental arrangement used a 3-mm-diameter collimated beam of 
800-nm, 150-fs, 1-kHz repetition rate light pulses from a Tirsapphire laser and amplifier system to 
illuminate one side of the slab sample. An ultrafast gated intensified camera system that provides a 
minimum FWHM gate width of 80 ps recorded the 2-D intensity patterns of the light transmitted through 
the opposite side of the slab. The gate position was varied in steps of 100 ps over a 5-ns range to obtain 
a sequence of 2-D transmitted light intensity patterns of both less-scattered and multiple-scattered light 
for image reconstruction. The inversion algorithm is based on the diffusion approximation of the 
radiative transfer theory for photon transport in a turbid medium. It uses a Green’s function pertur¬ 
bative approach under the Rytov approximation and combines a 2-D matrix inversion with a one¬ 
dimensional Fourier-transform inversion to achieve speedy 3-D image reconstruction. In addition to the 
lateral position, the method provides information about the axial position of the object as well, whereas 
the 2-D reconstruction methods yield only lateral position. © 1999 Optical Society of America 
OCIS codes: 100.3190, 170.3010, 170.3880, 170.6920, 290.7050. 


1. Introduction 

The past decade has witnessed a rapid growth of 
research activities in optical imaging of objects lo¬ 
cated inside turbid media. 1-12 This interest derives 
from potential practical applications of optical imag¬ 
ing in such diverse areas as detection of tumors in the 
human body, monitoring of aerosols and hydromete¬ 
ors in the atmosphere, and tracking of fuel droplets in 
the front nozzle of jet engines. Among these areas, 
biomedical applications have attracted the most at¬ 
tention because safe, noninvasive, and affordable im- 
aging modalities with diagnostic ability can be 
developed using near-infrared light. 2 An intensely 
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pursued area is optical mammography in which near- 
infrared light is used to screen for breast cancer, a 
major and increasing health problem that affects one 
out of every eight women in the United States. 13 
X-ray mammography, the most commonly used 
breast imaging modality, is not suitable for imaging 
young dense breasts, may not distinguish between 
malignant and benign tumors, and frequent exposure 
to ionizing x rays could be harmful if used for routine 
screening. The impetus for developing optical mam¬ 
mography is to overcome some of these limitations of 
x-ray mammography and to achieve diagnostic abil¬ 
ity by exploiting potential spectroscopic differenc¬ 
es 14 * 15 between the tumor and the surrounding 
tissues. 

In its simplest form, optical imaging involves illumi¬ 
nating an object with light of appropriate wavelength 
and looking for its shadow in the transmitted light. 
Some difference in optical properties between the ob¬ 
ject and the surrounding medium is necessary for the 
formation of a shadow image. However, strong scat¬ 
tering of light by biological tissues severely degrades 
the image quality, and for sufficiently thick tissues, 
completely buries it in the background noise. Time- 


1 July 1999 / Vol. 38, No. 19 / APPLIED OPTICS 4237 


22 




resolved, 3 ’ 5 ’ 61112 and continuous-wave (cw) 7 8 methods 
have been developed to sort out image-bearing ballistic 
and snake photons from the image-blurring multiple- 
scattered photons. 16 These methods improve the im¬ 
age quality to different extents, but suffer from low 
signal-to-noise ratio when used for imaging through 
thick, highly scattering samples. Inverse image re¬ 
construction (HR) methods make use of the measured 
scattered light intensity around the object, knowledge 
(or estimate) of optical properties of the turbid medium 
and the object inside it, and mathematical models to 
construct images of the object. 4 ’ 91217 " 24 An HR 
method requires extensive measurement of light in¬ 
tensities around the object and a long computation 
time for obtaining an image. Continuous-wave or 
frequency-modulated light, as well as ultrashort light, 
pulses may be used to accumulate data for HR. 

In this paper we report on a novel HR approach that 
makes use of experimentally measured time-sliced 
two-dimensional (2-D) transmitted light intensity dis¬ 
tribution, I(x, y , tg), to reconstruct three-dimensional 
(3-D) images. Time-sliced imaging is an extension of 
the idea of time-resolved early-light imaging except 
that picoseconds-duration temporal slices of the trans¬ 
mitted light covering both early-arriving less-scattered 
and late-arriving multiple-scattered photons are used 
to obtain a sequence of 2-D intensity patterns I(x, y , t t ). 
Transmission measurements are carried out using a 
gated image intensifier that provides an approxi¬ 
mately 80-ps duration time gate whose temporal posi¬ 
tion can be varied over a 20-ns range to collect a 
sequence of 2-D transmitted light intensity distribu¬ 
tion I(x, y, ti). The inversion algorithm uses a 2-D 
matrix inversion with a one-dimensional (1-D) 
Fourier-transform inversion based on symmetry in cy¬ 
lindrical coordinates that greatly reduces computation 
time compared with direct 3-D inversion. 

Our approach is a significant advance in optical to¬ 
mographic imaging of turbid media. First, use of 
time-resolved data enables reconstruction of 3-D im¬ 
ages that provide considerable information about 
object location in three dimensions. Depth informa¬ 
tion may not be available from direct transillumina¬ 
tion methods, or even x-ray mammography that yields 
only 2-D images. It has been demonstrated recently 
that HR approaches using cw measurements may not 
be adequate for a unique reconstruction of the internal 
optical characteristics of a scattering medium even 
when multiple sources and detectors are used. 25 
Shifting of the positions of the source detector pairs 
may lead to 3-D image information, but it would make 
data acquisition highly time-consuming. Second, 
time-sliced data provide a wealth of information over a 
broad frequency range that is essential for 3-D image 
reconstruction as compared with frequency-domain 
schemes that commonly use several discrete frequen¬ 
cies. Frequency-domain methods may, in principle, 
scan over a broad continuous range of frequencies to 
generate data that are equivalent to time-domain 
methods. However, accumulation and processing of 
data over such a broad frequency range (several giga¬ 
hertz) would be extremely time-consuming. More- 
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over, currently available sources cannot provide 
adequate power at such high frequencies. Use of 
time-resolved measurements is an effective way to col¬ 
lect data with adequate information content for 3-D 
image reconstruction. 18 Another advantage of the 
present scheme that makes use of time-sliced 2-D light 
intensity patterns over earlier time-domain HR 
schemes 17 ’ 18 is speedy data acquisition. Those earlier 
schemes involving measurement of time-resolved pro¬ 
files using multiple fibers and single-element detectors 
were time-consuming, expensive, and required compli¬ 
cated synchronization schemes. In the present 
scheme using a gated charge-coupled device (CCD) 
camera for detection, each pixel of the 384 X 288 pixels 
sensing element can be considered to be equivalent to 
a detector in the multidetector configuration men¬ 
tioned above, and a 2-D spatial distribution of trans¬ 
mitted light intensity is obtained in every frame 
recorded by the camera. Use of only the transmitted 
light measured with a single CCD camera makes data 
acquisition simpler and faster than data acquisition 
that requires measurements of light intensity around 
the object using multiple detectors. Finally, use of a 
2-D matrix inversion with a 1-D Fourier-transform 
inversion in the algorithm, together with fast time¬ 
gated data acquisition, enables reconstruction of 3-D 
images within a short time, a key requirement for 
real-time clinical applications. Although the linear 
perturbative approach has been used before, this use of 
a 2-D matrix inversion with a 1-D Fourier-transform 
inversion for speedy image reconstruction is a new 
salient feature of our approach. 

As a proof of concept, we used the method to recon¬ 
struct images of an absorbing object placed inside a 
tissue-simulating model turbid medium. We obtain 
good localization of the object in the lateral dimen¬ 
sions. Spatial resolution along the axial direction is 
promising, but needs improvement. The remainder 
of this paper is organized as follows. Section 2 de¬ 
tails the experimental arrangement, sample charac¬ 
teristics, and data-acquisition scheme. The HR 
algorithm is presented in Subsection 3.A. Because 
the inversion problem is ill-posed, some form of reg¬ 
ularization is needed to extract useful image infor¬ 
mation. The degree of regularization is typically 
controlled by parameters that are usually chosen em¬ 
pirically. The L-curve method 26 that we use to se¬ 
lect optimal regularization parameters is described in 
Subsection 3.B. The formalism is used to recon¬ 
struct images of an absorbing object inside a turbid 
medium and the results are presented in Section 4. 
The results of reconstruction using simulated data 
presented in Subsection 4.A is compared with that 
obtained using experimental data presented in Sub¬ 
section 4.B. Implications of these results for obtain¬ 
ing tomographic images of biological tissues and the 
scope of this HR approach are discussed in Section 5. 

2. Experimental Method and Materials 

The experimental arrangement for picosecond- 
duration time-sliced measurement of transmitted 
light intensity is displayed schematically in Fig. 1(a). 
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Fig. 1. (a) Schematic diagram of the experimental arrangement 

used for time-sliced transmitted intensity measurements, (b) 
Geometric arrangement used to describe the reconstruction algo¬ 
rithm showing the object and image planes and the object in an 
arbitrary position in cylindrical coordinates, (c) Object location in 
a lateral plane parallel to and 15 mm away from the exit plane. 
The 10-mm side of the sample was inclined at an angle of approx¬ 
imately 15° with the vertical (x axis). This is the object position 
that was used in the experiment. 


The detail of the cylindrical sample cell geometry is 
presented in Fig. 1(b). The location of the object on 
a plane parallel to and 15 mm toward the source from 
the exit plane (45 mm from the input plane) is shown 
in Fig. 1(c). The scattering medium was a suspen¬ 
sion of Intralipid-10% (Kabi Pharmacia Inc., Clayton, 
North Carolina) in water. Intralipid-10% is a fat 
emulsion used clinically as a nutrient and in research 
as a phantom to investigate light propagation in tis¬ 
sues. The concentration of Intralipid-10% suspen¬ 
sion was adjusted to provide an estimated reduced 
scattering coefficient jjl/ of 0.4 mm -1 and an absorp¬ 
tion coefficient jx a of 0.02 mm -1 at 800 nm. 27 The 
Intralipid suspension was held in a 60-mm-long and 
200-mm-diameter cylindrical Plexiglas cell. The 
cell thickness was chosen to be equivalent to that of 
an average breast under the extent of compression 
commonly used in x-ray mammography. 


The object to be imaged was a 3 X 3 X 10 mm 3 
rectangular parallelepiped made of aluminum and 
painted black. It was suspended on axis at a dis¬ 
tance of 15 mm (z = 45 mm) from the exit plane and 
absorbed most of the light incident on it. The cell 
was illuminated through the input plane with 150-fs 
duration, 1-kHz repetition rate pulses of800-nm light 
from a Ti;sapphire laser and amplifier system. 28 
The 3-mm-diameter collimated laser beam with an 
average power of approximately 100 mW was inci¬ 
dent along the axis (z axis) onto one of the flat faces 
(input plane) of the cylindrical cell. 

The 2-D intensity distribution of light emergent 
from the opposite endface of the cell (exit plane) was 
collected by a 24-mm focal-length ff 2.8 camera lens 
with a 84° angle of view and was recorded by the 
ultrafast gated intensified camera system (UGICS). 
The UGICS is a compact gated image intensifier unit 
that is fiber optically coupled to a CCD camera. It 
provides an electronic time gate whose FWHM dura¬ 
tion could be adjusted to a minimum of approxi¬ 
mately 80 ps. The gate position could be varied over 
a 20-ns range with a minimum step size of 25 ps. 
The signal recorded by the system at a particular gate 
position is a convolution of the transmitted light 
pulse with the gate pulse centered on the gate posi¬ 
tion. The collection optics was adjusted to capture 
light transmitted through the exit plane, and no at¬ 
tempt was made to use the scattered light emergent 
through the curved surface of the cell. The data 
then consisted of a sequence of time-resolved 2-D 
intensity distribution integrated over the gate dura¬ 
tion at different gate positions, I(x, y, t t ). Because 
I(x, y, t t ) is the 2-D spatial intensity distribution 
integrated over a 80-ps slice of time for the gate po¬ 
sition t if we refer to it as a time-sliced intensity pat¬ 
tern (or intensity distribution). These time-sliced 
intensity patterns were recorded by the CCD camera 
and displayed on a personal computer. 

In the experiment reported here, the gate width 
was adjusted to 80 ps and I(x, y , t t ) were recorded 
over a 5-ns range by varying the gate pulse position in 
steps of 100 ps. The zero time was taken to be the 
instant when the incident light pulse entered the 
front surface of the cell. I(x, y, t t ) provided the ex¬ 
perimental data for HR. 

Experimentally recorded typical 2-D intensity pat¬ 
terns of transmitted light with the object inside the 
medium are shown in Figs. 2(a) and 2(b) for gate 
positions of 700 ps and 1500 ps, respectively. The 
1500-ps pattern is spatially more spread out than the 
700-ps pattern because the electronic time gate col¬ 
lects light that is more diffuse and hence spatially 
spread out at the 1500-ps position than that at the 
700-ps position. Intensity patterns recorded with 
later gate positions spread out even further. Be¬ 
cause the temporal profile of the transmitted light 
intensity peaks around 2 ns, the intensity at the cen¬ 
ter of the 1500-ps pattern is higher than that of the 
700-ps pattern. No shadow image of the object was 
apparent in either of the intensity patterns. All the 
intensity patterns measured over the 5-ns range were 
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(a) t = 700 ps 


(b) t= 1,500 ps 



Fig. 2. Experimentally measured 2-D intensity distribution of 
light transmitted through the scattering medium with the object 
for a gate position of (a) 700 ps and (b) 1500 ps. (c) The object was 
located approximately on the axis of the cylindrical cell at a dis¬ 
tance of 15 mm from the exit plane. 


transferred to a Silicon Graphics Origin-2000 com¬ 
puter for use as experimental data for 3-D HR. The 
computer has four parallel 195-MHz CPU’s, but only 
one was used in this research. 

3. Reconstruction Formalism 

A. Algorithm 

The time-sliced HR formalism presented here is 
based on the diffusion approximation of the radiative 
transport theory for photon migration in a scattering 
medium. 9 The diffusion equation in a turbid me¬ 
dium is given by 

+ |x a (r)c - V[D(r)cV]j/(r, t) = S(r, t), (1) 

where I is the photon density (photons/cm 3 ), <S is the 
source strength [photons/(cm 3 -s)], |x 0 is the absorp¬ 
tion coefficient, D = l/{3[|x a + (1 - g) p. s ]} is the 
diffusion coefficient, is the scattering coefficient, g 
is the scattering anisotropy factor, and c is the speed 
of light in the medium. We are interested in the 
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following problem: Given the time-sliced intensity 
pattern transmitted through the media with the ob¬ 
ject I{x, y, t t ), and without the object I 0 (x, y, £,)> we 
want to reconstruct a map of the change in optical 
parameters, the absorption coefficient in this case. 
The medium without the object is referred to as the 
reference medium. 

The formulation of the forward problem uses a slab 
geometry (0, z 0 ) formed by cylindrical coordinates (p, 
<p, z), as shown schematically in Fig. 1(b). A point 
source (or a cp-symmetric finite source) is assumed to 
be located (or is centered) at (0, 0, 0). The detector 
system is located on the z 0 plane and is mapped onto 
the CCD camera by a lens. The p boundary is as¬ 
sumed to be at infinity. However, it can be replaced 
by a finite boundary using a more complex Green’s 
function. We use a Green’s function perturbative 
approach under the Rytov approximation 4 - 29 to pro¬ 
duce a linear inverse algorithm. 9 - 17 The Rytov ap¬ 
proximation is used because it has been shown to 
provide a more accurate reconstruction of the absorp¬ 
tive properties than the Bom approximation. 29 The 
reference medium, i.e., the medium without the ob¬ 
ject, is assumed to be uniform. Our formulation re¬ 
quires only the surrounding medium and not the 
object to be cylindrically symmetric, a condition that 
can be arranged in an experiment. For example, in 
the case of breast imaging the breast may be placed in 
a cylindrical cell containing an index-matching liquid 
with optical properties that are similar to the average 
values of optical properties of a real breast. The aim 
of the reconstruction process would then be to obtain 
a spatial map of the difference in optical properties 
from those average values. Thus the assumption of 
cylindrical symmetry for the surrounding medium is 
not a serious limitation for the applicability of the 
formalism. 

The forward problem can be represented in matrix 
form as 

Y = WX, (2) 

where Y is a vector represented by a column matrix 
whose elements are measured changes in transmit¬ 
ted intensity, X is a vector whose elements are optical 
parameter (absorption or scattering coefficient) 
changes that are due to the presence of the object, and 
W is the weight function, 9 which, in the linear inver¬ 
sion case, is related only to the reference medium. 
The vector Y has M elements, determined from the 
experimental data as Y t = -log {IJI 0 i), where I t and 
I oi are the measured intensities with and without the 
object, respectively, and M = (n t X n r X n^,) where n t 
is the number of time slices and n r and n, j, are the 
numbers of radial and angular sections, respectively, 
that the detector plane is divided into. X has N 
elements, Xj = cA|x 0 (rj) or cAZXr,), where N is the 
number of voxels that the sample cell is divided into. 
W is an M X N matrix that, under the assumptions of 
a cylindrical boundary and a uniform reference me¬ 
dium, satisfies the (p-rotation invariance, i.e., it is a 
function of cp - <p d , where <p d is the angle coordinate of 
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a point on the detector plane and 9 is the angle co¬ 
ordinate of a voxel. It should be pointed out that the 
reference medium (not the object) needs to be cylin- 
drically symmetric. For time-resolved absorption 
tomography, W is written as 


W 3 -d(p d , P, 9 “ 9 d*) = 


Ay 


i 


dT 


G P °(Pd, t)G°(z 0 , 0, t ) 

1 


■ivDc(t - t) 


X exp 


-p 2 - Pd 2 + 2 PPd cos(ip - <p d ) 


4 Defy - t ) 

X G z \z , z 0 ,t - t)G>, t)G,°(z, 0, T). (3) 

In Eq. (3), G z °(z, 0, t) is a 1-D slab Green’s function for 
the reference medium with source at z = 0 and is 
written as 

G z °(z, 0, t) = (4t TDct)~ l/2 exp(-z 2 /4Dct) 

+ 2 (4irZ)c£)~ 1/2 exp[—(2: - Zi) 2 /4Dct\, 


where i denotes the fth image source that is due to the 
slab geometry. 

G p °(p, t) is a 2-D Green’s function with the source 
located (or centered) at p = 0. For a point source, it 
is of the form 


(4) 


G '” (p ' * i^Dct ‘“’’(joi) ’ 

and for a finite-size spot source, it is written as 

- 2 hj 

(5) 

where 7 0 is the zeroth-order modified Bessel function 
and f(p s ) is the intensity distribution of the incident 
light spot. 

After making a 1-D Fourier transform of W 3 _ D over 
9 ~~ 9d> we obtain .^-independent 2-D matrices, 
W 2 _ d (^) parameterized by k , with K the number of 
grid points in the Fourier k space. The indices of cp - 
cp rf in the 3-D real space correspond to indices of k in 
the Fourier space. W^^k) is an M' X N' matrix, 
with M' - (number of p d division) X (number of time 
slices), and N r ~ (number of p division) X (number of 
2 division). The normal form W 2 ^(k) T W 2 , D (k) is an 
N’ X N* matrix. We separately calculate the K in¬ 
verse matrices [W^ikyW^ik) + A(^)] -1 , k = 1, 
2 , ,K, where k(k) is a matrix for regularization. 

The computational complexity now is K times that of 
inverting a W 2 . D matrix, which is much less than that 
of inverting a W 3 . D matrix. 

Having experimental data Y(p rf , cp d , t ), we make a 
Fourier transform over cp d to obtain Y(p d , k , t ), k - 1, 
2, .. . , K. The technique for fast discrete Fourier 



Fig. 3. L curve used for choosing the regularization parameters 
for HR with simulated data. The object was located at a distance 
of 15 mm from the exit plane. 


convolution is used. We then use simple matrix 
multiplication to obtain the image in the k space: 

X(k) = Y(^) T W 2 . D (^)[W 2 . D (^) r W 2 . D (A) + A(k)]~\ 

k = 1, 2 ,... K. (6) 

After X(&), k = 1,2, ... K, is obtained, the inverse 1-D 
Fourier transform over k reproduces the X as a func¬ 
tion of 9 and gives the 3-D distribution of the changes 
in absorption (or scattering) coefficients X(p, 9, z). 

B. Regularization 

The inversion problem is ill-posed with a huge con¬ 
dition number, hence a regularization procedure is 
necessary to obtain a unique solution. We use a 
diagonal form of a regularization matrix A(k) = 
\(&)8y'. A convenient way to understand the effect of 
regularization and to choose the suitable regulariza¬ 
tion parameter is to use the L curve, 26 which shows 
the properties of the regularized solution X x that vary 
with X. Two properties of the X x are the norm ||XJ 
and the norm of the corresponding residual vector 
||WX x - Y(|. |fXj indicates the stability of the solu¬ 
tion and it decreases with an increase of X. ||WX x - 
Y|| indicates the error that is due to the introduction 
of regularization and it increases with an increase of 
X. As proved in Ref. 26, if ||X x || versus ||WX x - Y|| is 
plotted with different values of X, a curve similar in 
shape to the capital letter L is obtained. It suggests 
that a good regularization parameter is one that cor¬ 
responds to a regularized solution near the comer of 
the L curve because in this region there is a good 
compromise that keeps both ||X x || and ||WX x - Y|| rel¬ 
atively small. Figure 3 shows an L curve that we 
obtained for one of our simulation tests, which corre¬ 
sponded to W 2 _ D (k = 1). The simulation assumed a 
Gaussian-distributed noise level with a root-mean- 
square deviation of 5% (for brevity henceforth re¬ 
ferred to simply as 5% Gaussian noise) in each 
simulated value of Y. The range of the values of X 
varies from 10~ 8 to 10 8 . The values of X correspond¬ 
ing to positions near the comer of the L curve are 
approximately 0.5-1.0. Choice of a value of X be- 
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tween 0.5-1.0 in our inversion calculation led to a 
stable image of the object with the 5% Gaussian 
noise. 

4. Results 

A. Three-Dimensional Image from Simulated Data 

In the simulation, the thickness of the slab was taken 
to be z 0 = 60 mm. A point light source was assumed. 
The detector system was taken to be a circle with a 
radius R d — 30 mm. To choose grids of equal area, 
the angle 2tt was uniformly divided into 36 parts, 
whereas the radius R d was divided into 30 parts 
scaled such that the outer radius of the jth disk was 
Rj = {R/j/J) 1/2 , where j = 1, 2, 3 ,,J, with J = 
30. In each simulated temporal profile, intensities 
at 30 time slices uniformly distributed from 450 to 
1900 ps were taken. The number of sampling data 
was M — 30 X 36 X 30. Data both with and without 
the object were used to obtain Y t = — log(/;// G ), after 
adding the 5% Gaussian noise. The sample cell was 
divided into N — 10 X 36 X 20 voxels, with 20 divi¬ 
sions along the z direction and 10 divisions along the 
p direction in a 30-mm range in a way similar to that 
used for detector division. 

Elements of W were calculated using Eq. (3) with a 
background-reduced scattering coefficient |x s '° of 0.4 
mm -1 and an absorption coefficient \i a ° of 0.02 mm* 1 . 
Corresponding parameters of the object were p, s ' = 
0.4 mm* 1 and jx a = 0.4 mm* 1 . The regularization 
parameter X was taken to be 1.0 based on the L-curve 
analysis. Image of the absorbing object obtained us¬ 
ing Eq. (6) and the 1-D Fourier transform is shown in 
Fig. 3. The 20 consecutive circles in Fig. 4 represent 
the images at 3-mm intervals along the cylinder axis 
(z axis) from the input plane to the exit plane. A 
gray scale is used to display the spatial distribution of 
changes in the absorption coefficient, with black rep¬ 
resenting the maximum value. For simulation the 
object of 3-mm linear dimension was placed in a voxel 
located in the 15th division along the z axis (z = 45 
mm) that was approximately 15 mm away from the 
exit plane. 

The reconstructed image locates the object in the 
neighborhood of where it was placed, as can be seen 
from the concentration of black around the 15th cir¬ 
cle. The image [Ap, a (r) distribution] appears to be 
spread out both in the lateral and in the axial direc¬ 
tions, a consequence of the scattering and the diffu¬ 
sion approximation used for reconstruction. The 
FWHM of the image is approximately 6 mm in the 
lateral plane and approximately 12 mm in the axial 
direction. 

B. Three-Dimensional Image from Experimental Data 

In the HR using experimental data, we approximated 
the laser beam to be a point source located at the 
center of the z 0 — 0-mm plane. The p boundary of 
the container was assumed to be at infinity. The 
time-resolved 384 X 288 pixels image from the CCD 
camera, with and without the object, were trans¬ 
formed to provide 30 X 36 X 30 pieces of Y data, as 
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Fig. 4. Images reconstructed from simulated data when the ob¬ 
ject was located at a distance of 15 mm (z = 45 mm ) from the exit 
plane. The sequence of circles represents images at 3-mm inter¬ 
vals along the cylinder axis. The linear intensity scale spans the 
range 0-20. The maximum value of A|x a (r) is 0.04 mm -1 and the 
size of the object is 3 X 3 X 3 mm 3 


was done in simulation. The plot of ||XJ versus 
||WX x - Y|| for computing the regularization param¬ 
eter X showed only a small kink on the curve instead 
of a clear L shape. The norm of the residual vector 
||WX x - Y|| at kink position was approximately 2 or¬ 
ders of magnitude larger than that in the simulated 
case. This difference may be due to the deviation of 
the experimental conditions from the assumptions of 
the theoretical forward model or because of a higher 
noise level in the experiment than that used in the 
simulation. Based on our simulation results assum¬ 
ing higher noise levels, we believe that the difference 
between the simulated and the experimental results 
mainly comes from higher experimental noise. Pos¬ 
sible sources of experimental noise are pulse-to-pulse 
energy fluctuation in the laser beam and timing jitter 
in the gate position. Because we arrange a small 
object at the distance far enough from the source and 
the detectors, use of the diffusion equation with the 
Rytov approximation seems reasonable. Figure 5 
shows a 3-D image of the object as a sequence of 
frames at 3-mm intervals along the axial direction 
reconstructed using experimental data. The lateral 
position of the object is reconstructed to be near the 
cylinder axis in the image map, as one would expect. 
However, the image is spread out with a half-width of 
approximately 20 mm along the axial direction. The 
absorption coefficient distribution appears to be 
peaked more toward the detector end than the actual 
object position. 
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Fig. 5. Images reconstructed from experimental data with the 
object located at a distance of 15 mm (z — 45 mm ) from the exit 
plane. The sequence of circles represents images at 3-mm inter¬ 
vals along the cylinder axis. 


One of the edges of the object was placed on the axis 
for experimental simplicity only and should not be 
considered a limitation of the present HR approach. 
Reconstruction using simulated data produced a sta¬ 
ble solution even when the object ‘was located closer 
to the cylinder boundary than to the cylinder axis. 

5. Discussion 

Results of this study indicate that the time-sliced 
transmission measurements along with the diffusion 
tomographic HR method can provide a 3-D map of the 
location of the object inside a biological tissue- 
simulating turbid medium. A salient feature of the 
approach is how short a time it takes for inverse 
reconstruction. Our estimate shows that a 3-D im¬ 
age of 10 X 36 X 20 voxels can be obtained with a 
running time of less than 1 min on a Silicon Graphics 
Origin-2000 computer when only one of its four par¬ 
allel 195-MHz CPU’s is used in the computation. 
The speed is comparable with that of a 200-MHz 
Pentium CPU. Combined with a fast data- 
acquisition scheme that we developed, this algorithm 
shows the potential for near-real-time HR. 

The reconstructed image of the object has the lat¬ 
eral position that one would expect from the experi¬ 
mental arrangement. The axial position is less 
certain, as the change in absorption coefficient A|x (t (r) 
is more spread out in the axial direction. Compari¬ 
son with image reconstructed from simulated data 
indicates that better axial localization is obtained 
using simulated data than that from the experiment. 
The position of the object in the reconstructed image 


gets shifted toward the detector plane for both the 
simulated and the experimental data. The shift is 
more pronounced for the image reconstructed using 
experimental data. 

To further investigate this apparent shift in the 
peak position and the extent of localization, we car¬ 
ried out inverse reconstruction with simulated data 
(a) for different axial positions of the object and (b) for 
different noise levels. Reconstructed images using 
simulated data with a 5% Gaussian noise level for 
object locations of 30 mm (z = 30 mm) and 45 mm 
(z = 15 mm) from the exit plane are shown in Figs. 
6(a) and 6(b), respectively. Figures 6(a) and 6(b) 
along with Fig. 4 indicate that, as the object is moved 
further away from the detector end, its localization in 
the reconstructed images becomes poorer. The ab¬ 
sorption coefficient change Ap, (r) spreads out both in 
the axial and in the lateral directions with a conse¬ 
quent decrease in its peak value as the object dis¬ 
tance from the exit plane increases. The result is 
expected because the further away the object is from 
the exit plane the higher is the probability for the 
image-bearing light to spread out and multiply scat¬ 
tered light to encroach into the region of a geometric 
shadow cast by the absorbing object. The overall 
effect is a smearing out of the shadow image. The 
diffusion approximation-based model used here is 
sensitive to that spreading and reconstructs images 
whose size and sharpness depend on the object loca¬ 
tion. 

The shift in the peak position of reconstructed 
A|x a (r) distribution toward the exit plane was ob¬ 
served for all three positions of the object. A quali¬ 
tative explanation for this shift can be given in terms 
of the dependence of imaging sensitivity on the object 
location. It is a common observation that for diff use 
light imaging, if the object is located on the exit plane 
it is imaged with minimal distortion and optimal 
sharpness. Image sharpness decreases and distor¬ 
tion increases as the object is moved further into the 
scattering medium. Imaging sensitivity is then a 
function that peaks on the exit plane and gradually 
decreases toward the input plane. The A|x a (r) dis¬ 
tribution in the reconstructed image would be pro¬ 
portional to a convolution of this function with A|x a (r) 
distribution of the object and would be shifted toward 
the exit plane. 

The spreading out of the Ap, a (r) distribution leads 
to a decrease in its amplitude in the reconstructed 
image as well because the magnitude of the inhomo¬ 
geneity gets distributed over a larger region of space. 
For a small object in three dimensions, the spreading 
of its image can dramatically decrease the amplitude 
of imaged optical parameters, in this case A|x a (r). 
This accounts for the peak amplitude of 0.04 mm -1 in 
the reconstructed image of Fig. 4, although the input 
value of the inhomogeneity was 0.4 mm " 1 
Another consequence of the imaging sensitivity be¬ 
ing peaked on the exit plane is the observed higher 
noise level in the z = 60-mm frame (exit plane) of the 
reconstructed images of Figs. 4, 6(a), and 6(b). The 
noise in the exit plane seems to be more pronounced 
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(b) 


Fig. 6. Images reconstructed from simulated data when the ob¬ 
ject was located at a distance of (a) 30 mm (z = 30 mm) and (b) 45 
mm (z = 15 mm) from the exit plane. The sequence of circles 
represents images at 3-mm intervals along the cylinder axis. 


as the object is moved toward the input plane, al¬ 
though the same 5% Gaussian noise level was used 
for all three cases. The apparent higher noise in Fig. 
6(b) is a consequence of the way the Ap, a (r) distribu¬ 
tion is displayed in these figures. The ratio of peak 
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value of A|ji a (r) to the magnitude of the noise level 
determines how prominently the noise will be dis¬ 
played in the scheme used in Figs. 4, 6(a), and 6(b). 
Because of the object location of 45 mm from the exit 
plane, the peak value of A|a a (r) is much smaller than 
that of the other two cases, the noise level appears to 
be more pronounced in the £ — 60-mm frame of Fig. 
6(b). Simulations with higher noise levels lend sup¬ 
port to this observation. With a high enough input 
noise level, the noise in the exit plane could be as 
pronounced as the peak value of A|x a (r) in the image 
of the object. 

The above observations of peak shift and exit plane 
noise with simulated data help explain the apparent 
lower localization and larger shift toward the detec¬ 
tor plane of the image from experimental data shown 
in Fig. 5. We ascribe these to the above-mentioned 
peaking of imaging sensitivity toward the exit plane 
and the higher noise level in the experiment com¬ 
pared with the 5% used in simulation. Sources of 
noise could be energy fluctuation of pulses from the 
laser amplifier, long-term drift of the laser output, 
jitters in the position and width of the time gate, and 
the CCD dark noise. Ample room exists for im¬ 
provement of the experimental noise and hence object 
localization. In this proof-of-the-principle experi¬ 
ment, our aim was to test if the algorithm provided 
meaningful 3-D images from time-sliced 2-D experi¬ 
mental data. 

One way of reducing the magnitude of peak shift 
and improving the localization of the object in the 
reconstructed image would be to use two sets of data 
for reconstruction—one taken with the source and 
detector positioned as described in the text and the 
second with the 2 : positions of the source and detector 
reversed. To test this concept, we carried out the 
reconstruction using the two sets of simulated data as 
mentioned above for the object position of Fig. 4, that 
is, 15 mm toward the source from the exit plane. 
The result, displayed in Fig. 7, shows improvement in 
object localization when compared with the image in 
Fig. 4. The improvement will be more substantial 
for an object located deeper into the scattering me¬ 
dium as a comparison of Fig. 7 with that of Fig. 6(b) 
reveals. 

The HR scheme presented above would be applica¬ 
ble to imaging of breast in vivo. The optical proper¬ 
ties, transport parameters, and thickness of the 
scattering medium used in this experiment are not 
too far removed from that of breast tissues. Al¬ 
though measurement on the breast with the lesion 
would provide the signal 7, the reference 7 0 can be 
obtained by using a uniform phantom with optical 
parameters that are equivalent to the average value 
of the parameters for the breasts. Another interest¬ 
ing possibility is to use measurements using the light 
of two different but nearby wavelengths. If a wave¬ 
length that is resonant with the lesion can be deter¬ 
mined, then measurements at that wavelength can 
be used as the signal and measurements with a non¬ 
resonant wavelength can be used as the reference. 
It should be noted that such lesion-specific wave- 
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Fig. 7. Images reconstructed from simulated data when the ob¬ 
ject was located at a distance of 15 mm (z — 45 mm) from the exit 
plane. Image reconstruction was carried out with data collected 
using the direct arrangement of both the source and the detector, 
as well as reversing their positions in the z direction. 


lengths are not well established yet. However, our 
spectroscopic imaging measurements indicate that if 
the imaging light is tuned into and out of the optical 
resonance of adipose breast tissues, then the adipose 
tissues could be highlighted in the transillumination 
image. 30 Combined spectroscopic and time-sliced 
measurements may provide much more information 
for tomographic image reconstruction than that 
available from single-wavelength x-ray mammogra¬ 
phy. Identification of fingerprint wavelengths for 
cancer diagnosis is an actively pursued area of cur¬ 
rent research. Imaging with such lesion-specific 
light, if identified, holds the promise for simultaneous 
detection, localization, and diagnosis of lesions, a 
unique potential advantage of optical mammogra¬ 
phy. 

One criticism sometimes raised about use of a ref¬ 
erence system is that measurements without the ob¬ 
ject in the reconstruction algorithm need to be 
addressed. Because our inverse algorithm is linear, 
the difference in the absorption coefficient between 
the object and the medium is relevant and the choice 
of a reference is arbitrary. In principle, it should be 
possible to reconstruct images from a set of time- 
sliced measurements on the object in the medium 
alone and take for the reference medium a uniform 
background characterized by the average values of 
optical properties. For example, in an optical mam¬ 
mography application, a uniform medium with val¬ 
ues of the optical properties similar to the respective 
average values of optical properties in the real breast 


could be considered a reference medium. The algo¬ 
rithm would try to estimate the difference between 
that assumed background distribution and measure¬ 
ment on the system with the object. The estimation 
can be carried out as long as the perturbation ap¬ 
proach remains valid. Alternatively, only measure¬ 
ments on the object within the medium using two 
different but nearby wavelengths of light could be 
used for the signal and reference, as discussed above. 

Our HR method presents a linear approach based 
on the Rytov approximation, and we were particu¬ 
larly interested in small objects, i.e., weak perturba¬ 
tion. One of the practical reasons for choosing small 
objects is that breast tumors iieed to be detected at an 
early stage of development when they are small in 
size and are more treatable. Also, the difference in 
optical properties between normal and cancerous tis¬ 
sues are generally not too large. The magnitude of 
perturbation depends on the size of both the object 
and the amplitude of changes in the optical proper¬ 
ties between the object and the surrounding medium. 
The combination is expected to be in the domain of 
small perturbation for early breast tumors. We con¬ 
sidered small objects with high optical contrast that 
seemed to satisfy the requirements for the linear ap¬ 
proach to hold as evident from the stable inverse 
solutions that we obtained. Linear approaches us¬ 
ing frequency-domain measurements were shown to 
yield inverse images for extended objects by other 
groups. 4 * 20 29 The present approach is expected to 
yield qualitatively adequate results for multiple ob¬ 
jects, and extended objects in the strong perturbation 
limit as well, although the details of the object shape 
and amplitude of optical parameter changes may not 
be quantitatively precise. However, at the current 
state of development, similar limitations are shared 
by other linear and nonlinear approaches as well. 
The advantages of the present linear approach are 
speed of reconstruction, ability to provide 3-D images, 
and generation of stable solutions without the need 
for artificial termination of iteration to avoid diver¬ 
gence that is needed for nonlinear methods. Unlike 
nonlinear iterative methods, 22 * 23 the linear approach 
is more susceptible to producing erroneous results for 
large and strong objects. 

In summary, even with the above-mentioned cave¬ 
ats and present limitations, the wavelength- 
dependent, time-sliced transmission measurements 
and the diffusion tomographic inverse algorithm 
show promise for providing 3-D images of objects in a 
turbid medium, such as a tumor in a breast within a 
clinically acceptable time frame. 
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ABSTRACT 

An inverse image reconstruction approach that makes use of an algorithm based on the diffusion approximation 
of the radiative transport theory and a sequence of picosecond-duration slices of transmitted two-dimensional (2D) 
light intensity distribution for fast 3D image reconstruction is presented. The results of simulation and experiment 
for a cylindrical geometry are presented. Object localization in the lateral dimensions is better than that in the 
axial direction. The observed difference in axial and lateral resolutions are analyzed by introducing the concept of 
“longitudinal sensitivity”. 

Keywords: Inverse reconstruction, time-sliced imaging, transillumination 

1. INTRODUCTION 

In the article, we present an inverse image reconstruction (IIR) approach that makes use of a sequence of picosecond- 
duration slices of transmitted two-dimensional (2D) light intensity distribution measured at one plane, instead of 
surrounding the medium, for fast 3D image reconstruction. This approach may allow for simpler 3D imaging of 
breast. 


2. OPTICAL IMAGING SYSTEM 

The experimental system, described in detail elsewhere, 1 makes use of 800 nm, 150/s, 1 kHz repetition-rate pulses 
from a Ti:sapphire laser for sample illumination and an electronic gated image intensifier coupled to a CCD camera 
to record time-sliced images. The scattering medium was a suspension of Intralipid-10% (Kabi Pharmacia Inc., 
Clayton, North Carolina) in water, adjusted to provide an estimated reduced scattering coefficient p! s = 0.4 mm -1 
and an absorption coefficient p a = 0.02 mm -3 at wavelength 800 nm. The Intralipid suspension was held in a 60 mm 
long and 200 mm diameter cylindrical Plexiglas cell. The object was a 3 x 3 x 10mm 3 rectangular parallelepiped 
made of aluminum and painted black. It was suspended on axis at a distance of 15 mm (z = 45 mm) from the 
exit plane. The gate width was adjusted to 80 ps and time-sliced transmitted intensity distribution, /(x,y,£ t ), were 
recorded over a 5 ns range by varying the gate position in steps of lOOps. 

3. ALGORITHM 

The inverse image reconstruction algorithm is based on the diffusion approximation of the radiative transport theory 
for photon migration in a scattering medium. For the linear inversion case, the forward problem can be written in 
matrix form: 

Y(Pd,<l>d,Zd,t) = Y, w (Pd,<j>d-<l>,Zd,t\p,0,z)X(p,<l>,z) (1) 

where V = - log(/// 0 ), I and / 0 are intensities with and without the inhomogeneity respectively, X = Ap a c or ADc 
depends on the type of the inhomogeneity, and W is the weight function. The convolution over <p in Eq. (1) can be 
eliminated by a Fourier transform in the 4> direction: 


2 


YkiPd, Zd,t) = w k(Pd,Zd,t\p,z)X k {p, z) (2) 

P.s 
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Figure 1 . (a) The cylindrical slab is uniformly divided into 20 slices along the z-axis. (b) Images reconstructed 
from simulated data when the object was located at z = 45 mm. (c) Images reconstructed from experimental data 
with the object located at z = 45 mm. The sequence of circles represents images at 3 mm intervals along the cylinder 
axis. 


where Y k , W k and X k are Fourier transforms of Y, W and X over <p respectively. The inverse problem is solved using 
a Tikhonov regularization method 2 with a diagonal form of regularization matrix A(Jfc). The suitable regularization 
parameters \{k) are chosen by the L-curve criterion. 3 In our calculation, the sample cell is divided into N r x N$ xN z = 
10 x 36 x 20 voxels. One inversion can be completed in 1 minute in a Silicon Graphics Origin200 workstation when 
the weight matrix is pre-computed. 


4. RESULTS 

4.1. Three-dimensional image from simulated data 

In the simulation, we assumed a cylindrical slab of thickness z 0 = 60 mm, and a disk detector of radius R d = 30 mm. 
The thickness of the cylindrical slab and its angle 27r was uniformly sliced into 20 and 36 parts respectively, while 
its radius R d was divided into 10 parts scaled such that each voxel has equal-volume (Figure 1(a)). Sampled data 
both with and without hidden objects were used to obtain Y = - log 7// 0 after added Gaussian distributed noise 
of ±5%. The background was assumed to have a reduced scattering coefficient p° s ' = 0.4 mm -1 , and an absorption 
coefficient p° - 0.02 mm -1 . The optical parameters of the inhomogeneity were p' s = pP s ' and p a = 0.4 mm -1 , which 
was put in a voxel at the 15th division along z-axis. 

The reconstructed image of Ap a is shown in Figure 1(b). The 20 consecutive circles in the figure represent the 
images at different position along the z-axis from input plane (z = 0) to exit plane (z = 60 mm). A gray scale is used 
to denote the relative value of the change of the absorption coefficient. The reconstructed image locates the object in 
the neighborhood of where it was placed. The object appears to spread out both in the lateral and axial directions, 
a consequence of the scattering and the diffusion approximation used for reconstruction. The image spreads out to 
approximately 6-7 mm in the lateral, (p, <p) plane, and to approximately 10-15 mm along the axial, z-direction. 

4.2. Three-dimensional image from experimental data 

In the inverse image reconstruction using experimental data, we approximated the laser beam to be a point source 
located at the center of the z = 0 plane. The p boundary of the container was assumed to be infinity. 

Figure 1(c) shows a 3D image of the object as a sequence of frames at 3 mm intervals along the axial direction from 
input plane to exit plane reconstructed using experimental data. The lateral position of the object is reconstructed 
to be near the cylinder axis in the image map, as the actual object was placed. However, the image is spread out over 
approximately 20 mm along the axial direction and peaked more towards the detector end than the actual object 
position. The worse result from experimental data indicates higher noise in experiment than 5% and existence of 
deviation between theoretical model and experiment. 


3 
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5. LONGITUDINAL SENSITIVITY ANALYSIS 

To understand why the axial resolution is poorer than lateral resolution for both simulated and experimental data, 
a one-dimensional model for an infinite uniform highly scattering medium with a diffusion coefficient D and an 
absorption coefficient ix a was studied. A source of plane-wave pulse with intensity S5(z + d)6(t) was located at 
z = -d < 0, and a detector at z = d. An inhomogeneity in the form of an absorptive infinite layer of infinitesimal 
thickness A 2 with diffusion coefficient D and absorption coefficient jj. a + Ap a was placed at 2 where -00 < 2 < 00 . 
The time-resolved fluxes received by the detector, without and with the inhomogeneity, are I 0 {t,z) and I\(t,z), 
respectively. The change of flux for an absorptive inhomogeneity is given by: 4 


A/(f, 2 ) = I 0 (t,z) - h{t,z) = SAfi a Az 



drG(d,z,t — t)G(z, —d,r) 


} 


(3) 


where G(zi,z 2 ,£) is the one-dimensional Green’s function. 

We define a sensitivity function as the term inside {} of Eq. (3) which is the change of flux measured by the 
detector at time t when an absorption layer of unit thickness and of unit change of absorption coefficient is placed 
at position z for an incident pulse of unit intensity. The sensitivity function is found to be 5 

Wmf(<, z) = e — erfc(amax(| 2 |,d)) (4) 

where a~ r = y/Dct is the characteristic length. W' m r(t,z) is totally independent of z when the absorptive layer is 
between the detector and the source. Therefore, for one-dimensional case, reconstruction using only transmission 
measurements is not sensitive to axial position of the inhomogeneity. In a 3D transmission case, such as that in 
section 3, it is not so severe but limits the longitudinal resolution, as we have observed in section 4. 

6. CONCLUSION 

An inverse image reconstruction approach that uses a sequence of time slices of two-dimensional light intensity 
distribution, and a fast algorithm based on the diffusion approximation of the radiative transport theory is presented. 
Results of reconstruction with simulated and experimental data are also presented. The image spreads to about 
6-7 mm in the lateral, and to approximately 10-15 mm along the axial direction for one voxel inhomogeneity in 
simulation while the image spreads out to approximately 20 mm along the axial direction for experimental data. The 
difference in lateral and axial resolutions may be explained by the concept of “longitudinal sensitivity”. 
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Fig. A schematic diagram of the time-sliced imaging arrangement with an optical Kerr gate using a Ti: sapphire laser 
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Fig. A schematic diagram of the time-sliced imaging arrangement with an optical Kerr gate using the Crrforsterite laser 
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Appendix 7 



Fig. A schematic diagram of the experimental arrangement for spectroscopic imaging. 
(A = aperture; P = polarizer; S = sample) 
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